Survival of the Aligned: Ordering of the Plant Cortical Microtubule Array 
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The cortical array is a structure consisting of highly aligned microtubules which plays a crucial role 
in the characteristic uniaxial expansion of all growing plant cells. Recent experiments have shown 
polymerization-driven collisions between the membrane-bound cortical microtubules, suggesting a 
possible mechanism for their alignment. We present both a coarse-grained theoretical model and 
stochastic particle-based simulations of this mechanism, and compare the results from these comple- 
mentary approaches. Our results indicate that collisions that induce depolymerization are sufficient 
to generate the alignment of microtubules in the cortical array. 

PACS numbers: 87.16.Ka, 87.16.ad, 87.16.af, 87.16.Ln 



Microtubules are a ubiquitous component of the cy- 
toskeleton of eukaryotic cells. These dynamic filamen- 
tous protein aggregates, in association with a host of mi- 
crotubule associated proteins (MAPs), are able to self- 
organize into dynamic, spatially extended stable struc- 
tures on the scale of the cell I ll . In contrast to the more 




Figure 1: Transverse cortical array in an etiolated dark-grown 
Arabidopsis thaliana hypocotyl cell with fluorescently labeled 
microtubules. Image courtesy of Jelmer Lindeboom, Wagenin- 
gen University. 



commonly studied animal cells, plant cells are encased in 
a cellulosic cell wall, and generally only expand along a 
single well-defined growth axis. A crucial component in 
this anisotropic growth process is a plant-unique micro- 
tubule structure called the cortical array [2] . This struc- 
ture consists of highly aligned microtubules attached to 
the inner side of the cell membrane and oriented trans- 
versely to the growth direction (see Fig. 1) and estab- 
lishes itself in a period of about one hour after cell di- 
vision. The cortical array has two particular features, 
both related to the fact that the microtubules are bound 
to the cell membrane 0, 01 : (i) it is effectively a 2- 
dimensional system and (ii) the cortical microtubules do 
not slide along the membrane, so the only displacements 
are caused by the ongoing polymerization and depolymer- 
ization processes intrinsic to microtubules. As a conse- 
quence of these two constraints, the so-called plus end of 
a growing cortical microtubule can 'collide' with another 
microtubule. Recent experiments [5| have shown that 
these collisions indeed occur and can have three possi- 
ble outcomes whose relative frequency is determined by 
the angle between the microtubules involved (see Fig. [2^) . 
The first option is that the incoming microtubule changes 
its direction and continues to grow alongside the micro- 
tubule it encountered, an outcome that is predominant 



at smaller angles and is known as 'zippering'. The second 
option is the so-called 'induced catastrophe', in which the 
incoming microtubule switches to the shrinking state. Fi- 
nally, there is a possibility that the incoming microtubule 
simply 'crosses over' the obstacle, continuing to grow in 
its original direction. 

In this Letter we address the question of whether, as 
has been posited by Dixit and Cyr ^ , these interactions 
are sufficient to explain the alignment of microtubules in 
the cortical array. To do so we construct a model for the 
microtubule dynamics and interactions, and evaluate it 
using two complementary approaches: a coarse-grained 
theory and particle-based simulations. The theory al- 
lows us to reduce the size of the model parameter space 
by identifying the relevant control parameter of the sys- 
tem and establishes the criteria for spontaneous symme- 
try breaking to occur. The simulations explicitly con- 
sider the stochastic dynamics of individual microtubules, 
and are thereby able to test the validity of the theory. 
The simulations can also be extended to include known 
other contributing effects such as minus-end treadmilling 
and microtubule severing proteins, but here we focus on 
a minimal version of the model that can be addressed 
using both the theoretical and simulation approaches in 
order to establish a reference system and test the general 
hypothesis of 0]. 

Our model differs from existing models for 2D orga- 
nization of filamentous proteins in two important ways. 
Firstly, in most of these models the filaments are both 
free to rotate and translate as a whole 0, S, 0, 0, U | , 
which is inconsistent with the experimental observations 
on the cortical array. Secondly, our model explicitly takes 
into account the dynamic instability of the individual mi- 
crotubules, providing the potential for intrinsic stabiliza- 
tion of the microtubule length distribution. This differs 
from the model by Baulin et al. [l^ in which determinis- 
tically elongating microtubules stop growing only while 
obstructed by other microtubules. The lack of an intrinsi- 
cally bounded length most likely precludes the existence 
of stable stationary states in their simulations. 

For the intrinsic microtubule dynamics in our model, 
we use the standard two-state dynamic instability model 
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y in which each microtubule plus end is assumed to be 
either growing with a speed or shrinking with a speed 
v~ . This plus end can switch stochastically from growing 
to shrinking (a so-called 'catastrophe') with rate rc, or 
from shrinking to growing (a so-called 'rescue') with rate 
in a process known as dynamic instability. New micro- 
tubules are nucleated isotropically and homogeneously 
with a constant rate rn- The microtubule minus ends are 
assumed to remain attached to their nucleation sites. 




Figure 2: a) schematic overview of the included effects and 
parameters in the model, b) relative frequency of collision 
outcomes as a function of angle of incidence used in our model. 



Because the persistence length Ip of microtubules is 
long [r^ mm) compared to the average length of a micro- 
tubule 10/im) and thermal motion is inhibited by the 
attachment to the plasma membrane, microtubules are 
modelled as straight rods with kinks at positions where 
a zippering event has occurred. A microtubule therefore 
consists of a series of connected segments to which we 
assign an index i, starting at i = 1 for the segment at- 
tached to the nucleation site. In light of the available 
evidence, we assume that the angle-dependent collision 
outcome probabilities (zippering), Pc (induced catas- 
trophe) and Px (crossover) are independent of the polar- 
ity of the microtubules and are therefore fully defined on 
the interval [0, ■!]. 



We first analyze this system using a coarse-grained the- 
ory, in which we consider densities of microtubule seg- 
ments instead of individual microtubules. From the out- 
set we assume the system is, and remains, spatially ho- 
mogeneous, and later restrict ourselves to steady-state 
solutions. Because microtubules are nucleated isotropi- 
cally and can change their orientation after each zipper- 
ing event, we introduce separate densities for each seg- 
ment index i. Furthermore, length changes and collisions 
can only occur in segments that contain the microtubule 
plus end. Therefore, we further distinguish the active 
segments, containing either a growing (-I-) or shrinking 
(-) plus end, and the inactive (0) segments that form the 
'body' and tail of the microtubule. Our variables are thus 
the areal number densities mf (/, 9, t) of segments in state 
a £ {0, — , +} with segment index i, having length I and 
orientation 9 (measured from an arbitrary axis) at time t. 
From these, we compute the total length density k{9, t) 
as 

k{9,t)^y2 dl I [m+{l,9,t) + {l,9,t) + m°{l,9,t)] . 

(1) 

The segment densities obey a set of evolution equations 
that can symbolically be written as 

dtmf{l,9,t) =$g„„[m+] -I- ^rescue[m~] - ^sp.cat.[rn:l] 

- ^ind.cat.iml, k] - ^^,p[mf,k] (2a) 

dtm~{l,e,t) =^ shrink ["^rl - ^rescneKi"] + ^ sp.cat.[mf] 
+ ^ind.cat.[m:l, k] + ^react.[m'l,mY^j^,k] 

(2b) 

dtm°{l,9,t) =(^^,p[m+,k] - $reQct.[TO+, m^^^, fc] (2c) 

The arguments in square brackets explicitly display the 
functional dependencies of the terms on the right hand 
side. Below, we explain each of these terms briefly, and 
refer the reader to [ist for a full derivation and an in- 
depth analysis. The dynamics of the active growing 
(-I-) and shrinking (— ) segments of microtubules unper- 
turbed by interactions are given by the standard sponta- 
neous catastrophe and rescue rates ^sp.cat.i'm^] — rcm^ 
and $rescue["^ ] = r-^m' , and the advective terms 
^gron,[m+] = ~'V+ ^ and $sft„„fe[m-] = v' ^ 
due to growth and shrinkage respectively Colli- 
sions between microtubules that lead to an induced 
catastrophe cause growing segments to switch to the 
shrinking state, at a rate given by ^ind.cat.[m'^ , k] ~ 
v+m+{9) J de' sinA9PciA9) k{9'), where M = \9 - 
6'\ is the collision angle and the geometrical factor 
sin A9 takes care of the collisional cross-section the den- 
sity of other microtubules present to the incoming one. 
Zippering events cause growing microtubule plus ends 
to change direction, converting previously growing seg- 
ments to the inactive state at a rate $2ip[m+. A:] = 
v+m+{9) J d9' sin A9Pz{A9) k{9'). Simultaneously, new 
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growing segments with an index i+1 are created, which is 
represented by the boundary condition m^-^(/ = 0,9) — 
J dl'^zip[mf , k]. This set of boundary conditions is com- 
pleted by a separate equation for i = 1, which represents 
the isotropic nucleation of new microtubules: v'^m'^^l = 
0,9) = rn/(27r). Finally, when a segment shrinks back 
to the point where it had undergone a zippering event 
in the past, a previously inactive segment can be reacti- 
vated into a shrinking state. Here we will not discuss the 
details of the rate ^reactimf , k], which contains a 

non-trivial history-dependence as a microtubule segment 
must "un-zipper" in the same direction the zippering seg- 
ment originally came from. We simply note that in the 
steady state Eq. (|^ requires that this rate is balanced 
by the zippering rate discussed above. 

In the steady state, the infinite set of equations ^ 
with the boundary conditions reduces to a set of four 
coupled non-linear integral equations. These relate the 
length density k{9) to the average segment length, active 
segment density and ratio between inactive and active 
segments, each being a function of the angle 9. It follows 
that, for given interaction probabilities Pc{9) and Pz{0), 
the remaining parameters can be absorbed into a single 
dimensionless control parameter G, defined as 



G = 



fn (w+ + V ) 



( H. - 



(3) 



Here we only consider the case G < 0, for which the 
length of the microtubules is intrinsically bounded even 
in the absence of collisions. In this case, the aver- 
age length of non-interacting microtubules is given by 
T — {rc/v+ — rr/v~)~^ and the control parameter G 
can be interpreted as G = —lo/l, implicitly defining an 
interaction length scale Zq. As G increases towards 0, the 
number of interactions between microtubules increases. 

For any value of G there exists an isotropic solution 
to ([2]), for which the total length density p — Jd9k{9) 
satisfies Iqp {cqIqp — 2G)^ = 8, where c„ denotes the n- 
th Fourier cosine coefficient of the product Pc{9) \sm9\. 
The isotropic length density is therefore an increasing 
function of the control parameter G that only depends 
on the induced catastrophes, and not on the probabil- 
ity of zippering. This can be understood by the fact 
that zippering only serves to reorient the microtubules, 
which has no net effect in the isotropic state. Although 
a stationary isotropic solution exists for all values of G, 
this solution is only stable for large negative values of G. 
As G increases, the number of interactions between mi- 
crotubules increases, until the isotropic solution becomes 
unstable. This happens at the bifurcation point G = G* , 
given by 



G* = (-2c2) 
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(4) 



We note that the location of the bifurcation point is deter- 
mined solely by the properties of the induced catastrophe 



probability Pc{9), and, like the density in the isotropic 
phase, does not depend on zippering. 
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Figure 3: Comparison between theoretical (solid lines) and 
simulation results (symbols) . The simulations were performed 
on a 80 pm x 80 p.m system with periodic boundary con- 
ditions. The spontaneous catastrophe rate was varied to 
probe different values of G: rc £ [4 x 10"^, 1.2 x 10"^] 
s^^. The nucleation rate was set to r„ = 0.003 /im^^s^^ 
and other parameters were taken from 4] (interphase BY- 
2 cells); v+ = 0.078 pm s-\ v~ = 0.164 pm s'^ , r,. = 6.8 
xlO"'' s~^. Measurements were performed after equilibrat- 
ing for 50,000s (a) or 250,000s (b); G was increased between 
measurements. The standard error of the mean is typically 
smaller than the symbols and is otherwise indicated by verti- 
cal bars. N =80(a), 40(b). 



To quantify the degree of alignment we use the stan- 
dard 2D nematic order parameter 5*2, defined as ^2 = 
ll"" Aee^'^k{9)\/ /q " dO k{9). The full bifurcation dia- 
gram can be computed by numerically tracing the or- 
dered solution branch from the bifurcation point, pro- 
vided that the products | sin6'|Pc(^) and | s\\\9\Pz{9) have 
finite Fourier expansions. We restrict ourselves to an ex- 
pansion up to cos 40. The coefficients are constrained by 
sin (O)Pc(O) = and sin (O)P^(O) = 0. In line with experi- 
mental observations Q we choose the remaining parame- 
ters such that Pc{9) is monotonically increasing to a max- 
imum at = 7r/2 and is maximally biased towards steep 
collision angles (see for other choices), and Pz{n/2) = 
0. The magnitudes of Pc{9) and Pz{9) is similar to that 
observed in experiments, and the crossover probability 
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is fixed by the requirement Pc{9) + Pz{0) + Px{0) = 1- 
The resulting interaction probabilities are illustrated in 
Fig. [2b- We argue that the apparent discrepancy with 
experiments, caused by setting Pz(0) = 0, is not very sig- 
nificant for the ordering transition, as collisions between 
near-parallel microtubules are infrequent and cause only 
slight changes of orientation in the case they lead to zip- 
pering. 

Given our choice for Pc{d), we have cq = 3/8 and 
C2 — -1/4 so that G* « —0.2. The results are rep- 
resentative for a large class of interaction probabilities 
with G* < 0. Higher modes do not affect the bifurca- 
tion point (U and appear to have only minor effects on 
the bifurcation diagram. Also, any changes to the over- 
all magnitude of Pc{0) and Pz{0) result only in a scahng 
of the G-axis. Comparing the computed solutions (solid 
lines) for systems with (Fig. and without ((2^) zip- 
pering, we note that zippering has only a minor effect 
on the ordering beyond the bifurcation point (see also 
[ist). This shows that the 'weeding out' of microtubules 
in the minority direction through induced catastrophes 
is by itself sufficient to explain microtubule alignment. 

In parallel with the coarse-grained theoretical ap- 
proach described above, we performed stochastic particle- 
based simulations of the interacting microtubules. Fig. [3] 
shows the resulting steady-state alignment as a function 
of G, for systems with and without zippering. In the 
simulations, the presence of zippering triggers the for- 
mation of microtubule bundles, in which aligned micro- 
tubules colocalize. In this case, we need to specify how 
the interaction probabilities Pc{S), Pz{0) and Px{0) de- 
pend on the number of microtubules that are present in 
both the incoming and encountered bundles. We investi- 
gate two extreme scenarios. In the first scenario (single 
collisions) a microtubule treats a collision with a bundle 
as a single collision, disregarding the other microtubules 
in both bundles. In the other scenario (multi-collisions) 
we implicitly construct an effective interaction by sam- 
pling from the distribution of all multiple collisions and 
their outcomes that can occur between an arbitrary mi- 
crotubule from an incoming bundle with the full set of 
microtubules in the target bundle (see Fig. [2k)- 

In the absence of zippering Fig. [3^, shows that the 
theoretical predictions and simulation results agree well. 
As expected, the agreement is less good when zippering 
is enabled (Fig. [3)3), because zippering leads to strong 
spatial correlations in the form of microtubule bundles, 
which are not accounted for in our mean-field-like theory. 
In the case of the 'multi-collision' interactions, the sim- 
ulations indicate a significantly larger tendency to align, 
whereas the system is less likely to align with 'single' in- 
teractions. However, in both cases the behavior remains 
qualitatively the same as the theoretical prediction and 
the alignment occurs over a similar range of G values. 

Finally we investigated the limit of weak interactions 
{Pc{0), Pz{0) ^ 1; data not shown) in which the discrep- 



ancies due to the mean-field nature of our model should 
decrease. Without zippering simulation results rapidly 
converge to the theoretical predictions. In the presence of 
zippering the results for the 'single' interactions deviate 
more strongly from the theory, because only a single col- 
lision is registered when a microtubule encounters a bun- 
dle, effectively decreasing the density of interactions. The 
'multi-collision' interaction however effectively accounts 
for the bundling, so that for progressively weaker inter- 
actions the transition between the isotropic and ordered 
states converges to the predicted bifurcation point. 

Our model of interacting cortical microtubules displays 
both isotropic and aligned phases and is based on ex- 
perimentally observed microscopic effects. The kinetic 
parameters appearing in the control parameter G may 
be regulated by the cell via MAPs, suggesting a mecha- 
nism for cellular control over creation, maintenance and 
suppression of microtubule alignment. Our results in- 
dicate that collision-induced microtubule catastrophes 
alone could establish alignment in the cortical array of 
plant cells. To what extent other known effects, such 
as microtubule treadmilling and severing, influence this 
mechanism, is a question we are currently addressing. 
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